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cotermination  of  the  constants  plus  one  more,  which  consists  of  a  single  rotation  about 
E.  It  is  noted  that  the  measurement  of  pure  thickness  vibrations  is  inadequate  for  the 
determination  of  the  temperature  derivatives  of  the  effective  piezoelectric  and  dielectric 
constants.  However,  the  solution  for  the  trapped  energy  resonator  is  presented  and 
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obtained  from  measurements  of  trapped  energy  resonance  in  thermally  compensated  cuts. 

'•’.'hen  the  measured  dac  >  is  available  all  coefficients  will  be  determined  from  the  analytical 
work  presented  here. 
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1 .  Introduction 


The  elastic,  piezoelectric  and  dielectric  constants  of  quartz  and 
their  behavior  with  temperature  have  a  significant  influence  on  the 
frequency  characteristics  of  various  quartz  devices  and  their  variation 
with  temperature.  Clearly,  the  more  accurately  these  quantities  are 
known,  the  more  precisely  quartz  device  structures  and  their  orientations 
can  be  described.  The  increased  precision  can  provide  advantages  in  both 
the  design  of  and  improvement  in  the  characteristics  of  resonant  quartz 
devices.  At  the  present  time  the  values  of  the  elastic,  piezoelectric  and 
dielectric  constants  are  those  of  Bechmann^  at  25^0  and  the  temperature 

derivatives  of  the  elastic  constants  that  are  used  are  those  of  Bechmann, 

") 

Ballato  and  Lukaszek".  The  accuracy  of  the  coefficients  currently  in  use 
has  been  subject  to  question,  in  particular  by  Kahan  ,  who  made  a 
statistical  comparison  of  the  best  measured  data  available  with  calcula¬ 
tions  based  on  different  sets  of  material  coefficients  and  concluded  that 
the  material  coefficients  and  their  temperature  derivatives  should  be 
redetermined. 

Before  proceeding,  it  is  important  to  note  that  since  the  existing 
temperature  derivatives  of  the  elastic  constants  of  quartz  were  evaluated 
from  the  data  using  the  linear  theory  of  elasticity,  which  can  only  be 
referred  to  the  temperature  dependent  intermediate  position  of  the  plate, 
they  are  the  temperature  derivatives  of  certain  effective  coefficients 
rather  than  the  fundamental  constants,  which  are  referred  to  a  fixed 
reference  position  at  one  temperature.  As  a  consequence  of  this,  the 

first  temperature  derivatives  of  the  fundamental  elastic  constants  of 

4 

quartz  were  subsequently  obtained  from  the  data  in  Ref. 2  within  the 


framework  of  the  proper  rotationally  invariant  nonlinear  thermoelastic 
description,  in  which  the  vibration  is  treated  as  a  small  linear  dynamic 
field  superposed  on  the  thermally  induced  static  biasing  state  and  referred 
to  the  temperature-independent  reference  position  of  the  plate.  Since  the 
vibration  is  referred  to  reference  coordinates,  the  mass  density  and  plate 
thickness  are  constants  independent  of  temperature  and  the  normal  to  the 
major  surfaces  of  the  plate  does  not  change  its  direction  with  respect  to 
the  principal  axes  of  the  quartz  crystal  with  temperature.  Since  in 
quartz  the  principal  axes  only  extend  and  contract  with  temperature  and 
all  others  skew,  the  actual  normal  to  the  intermediate  position  of  the 
surfaces  of  the  plate  changes  with  temperature.  This  change  in  normal 
(or  skewing  of  the  axes)  was  neglected  in  Ref. 2,  which  uses  the  linear 
description  referred  to  the  intermediate  coordinates.  This  is  the  primary 
reason  that  the  nonlinear  description,  which  permits  everything  to  be 
referred  to  the  unchanged  reference  coordinates,  has  a  significant 
advantage  over  the  commonly  used  linear  description.  However,  in  the 
determination  of  the  first  temperature  derivatives  of  the  fundamental 
elastic  constants  in  Ref. 4  only  the  elastic  (and  not  the  piezoelectric) 
solution  was  used  in  obtaining  the  temperature  derivatives  from  the  data 
and  the  temperature  derivatives  of  the  piezoelectric  and  dielectric 
constants  were  expressly  ignored.  Furthermore,  the  (rather  thick)  elec¬ 
trodes  on  the  quartz  plates  used  in  the  measurements  were  ignored  in  the 
treatment,  as  in  Ref. 2.  This  is  probably  the  major  source  of  inaccuracy 


in  the  existing  coefficients. 


The  vibration  solution  alluded  to  In  the  previous  discussion  is  the 

pure  thickness  solution”",  which  ignores  the  transverse  mode  shape  due  to 

the  finite  dimensions  of  the  electrodes  and/or  the  quartz  plate.  This 

procedure  is  reasonable  in  the  case  of  thermally  uncompensated  cuts  because 

the  iniluence  of  the  transverse  mode  shape  is  small  compared  to  that  of 

the  thickness  behavior.  However,  in  the  case  of  thermally  compensated 

cuts  such  as,  e.g.,  the  AT  and  SC  cuts,  which  are  the  most  important  in 

practice,  the  transverse  modal  behavior  is  of  crucial  importance  in 

determining  the  variation  in  frequency  with  temperature  because  for  those 

cuts  the  change  in  frequency  with  temperature  for  the  pure  thickness  mode 

of  interest  vanishes.  Furthermore,  the  temperature  dependence  of  the 

motional  capacitive  effect  of  the  driving  electrodes  on  the  quartz  plate, 

which  depends  on  the  temperature  derivatives  of  the  pertinent  piezoelectric 

and  dielectric  constants  for  the  thickness  mode  of  interest,  causes  the 

£ 

well-known  apparent  shift  in  angle  of  the  zero  temperature  cut  elcctroded 
quartz  plate.  In  a  calculation  of  the  temperature  dependence  of  the 
resonant  frequency  of  contoured  AT-cut  quartz  plates the  temperature 
dependence  of  the  motional  capacitive  effect  of  the  thickness  mode  of 
interest  had  to  be  estimated  from  temperature  measurements  on  AT-cut  quartz 
trapped  energy  resonators  with  large  electrodes  because  the  temperature 
derivatives  of  the  piezoelectric  and  dielectric  tensors  are  not  presently 
known.  In  addition,  because  of  the  inaccuracies  in  the  first  temperature 
derivatives  of  the  fundamental  elastic  constants  the  calculated  rotation 
angle  of  the  zero  temperature  AT-cut  unelectroded  flat  plate  is  -35°  15', 
which  is  referred  to  as  nominal  because  the  rotation  angle  of  the  actual 
cut  is  about  -35'”  21*  .  The  difference  of  6’  is  primarily  a  result  of  the 
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inaccuracy  caused  by  the  electrodes  that  were  on  the  plates  when  the 

measurements  were  made  and  were  not  considered  in  the  accompanying 
2  4 

analyses  ’  .  Similarly,  in  the  case  of  the  doubly-rotated  SC-cut  quartz 

g 

plate  using  the  same  6’  correction  to  the  9-angle,  we  have  found  with 

9 

the  aid  of  a  measurement  by  Warner  that  the  required  correction  to  the 
cp-angle  is  about  48',  which  is  quite  a  bit  larger  than  the  required 
correction  to  the  6-angle.  The  greater  error  in  the  co-angle  is  not 
surprising  because  it  relies  on  data  obtained  from  the  measurement^  of 
doubly-rotated  cuts.  Furthermore,  as  in  the  case  of  the  AT-cut,  the 
motional  capacitive  effect  due  to  the  driving  electrodes  on  the  SC-cut 
cannot  be  calculated  because  the  temperature  derivatives  of  the  piezo¬ 
electric  and  dielectric  constants  are  not  presently  known.  Again  the 

8  11 

effect  has  been  estimated  from  measurements  by  Lukaszek  on  that 
particular  cut.  However,  it  is  clearly  undesirable  to  perform  such 
measurements  on  every  zero  temperature  cut,  and  the  error  in  the  m-angle 
borders  on  the  intolerable. 

In  view  of  the  existing  situation,  the  elastic  piezoelectric  and 
dielectric  constants  of  quartz  and  their  first  four  temperature  deriva¬ 
tives  are  being  redetermined  at  a  fixed  reference  temperature  of  25°C. 

In  the  earlier  determination^  of  the  elastic,  piezoelectric  and  dielectric 
constants  of  quartz  only  thickness-excitation  of  thickness  vibrations  was 
employed  along  with  a  judicious  use  of  rods  and  contour  modes  of  plates 
and  the  then  existing  state  of  analytical  knowledge.  As  a  consequence, 
only  the  elastic  constants  c,,  and  c,,  could  be  determined  from  unrotated 

i  I  DO 

cuts  and  the  others  had  to  be  determined  from  some  singly-  and  some 
doubly-rotated  cuts  with  an  attendant  loss  in  accuracy.  Linder  this  program 
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both  thickness-excitation  and  lateral-excitation  of  thickness  vibrations 
are  being  employed.  This  should  result  in  a  significant  increase  in 
accuracy. 

m  order  to  determine  the  material  constants  from  measurements  of 

thickness  resonances  an  analysis  of  an  arbitrarily  oriented  quartz  plate 

driven  into  thickness  vibrations  by  either  thickness-excitation  or  lateral 

field  excitation  has  been  performed.  For  the  case  of  thickness-excitation 

the  analysis  is  restricted  to  orientations  for  which  the  three-coupled 

waves'  essentially  uncouple  and  one  dominates  the  vibration.  This  holds 

12  o 

for  all  orientations  except  those  tor  which  0  <  s  <  +  15  and  +22“  ^  <  30 

For  the  case  of  lateral  excitation  this  restriction  cn  orientation  holds 

even  chough  the  waves  are  not  coupled  because  they  are  almost  degenerate. 

The  analysis  reveals  that  the  four  constants  c,.,  c,.,  c,,  and  c,,  can  be 

determined  from  measurements  on  the  three  unrotated  cuts  and  that  the 

remaining  two  constants  c^  anc*  oan  determined  from  measurements 

on  rotated  Y-cuts.  Consequently,  no  doub ly- rota ced  orientations  are 

required  for  the  determination  of  the  constants.  Clearly,  this  should 

result  in  a  significant  increase  in  accuracy  both  because  of  the  simpler 

orientation  and  the  increased  directness  of  the  equations.  Since  the 

piezoelectric  constants  are  being  determined  from  the  measurement  of 

1 3 

successive  thickness  overcone  resonances  .  they  should  be  more  accurate 
than  in  the  earlier  work\  which  used  antiresonance  measurements,  the 
interpretation  of  which  requires  a  great  deal  of  insight.  The  dielectric 
constants  are  also  being  determined  from  the  same  measurements  of  over¬ 
tone  thickness  resonances. 


In  order  to  determine  the  temperature  derivatives  of  the  effective 

material  constants  a  perturbation  analysis  of  the  temperature  dependence 

of  the  resonant  frequencies  of  arbitrarily  oriented  quartz  plates  vibrating 

in  pure  thickness  modes  has  been  performed.  Since  as  already  noted  a 

14 

proper  rotationally  invariant  nonlinear  description  ,  which  enables  the 

equations  to  be  referred  to  a  fixed  reference  position  at  a  fixed  reference 

temperature  T  is  being  employed,  the  geometry  and  density  do  not  change. 

As  a  consequence,  the  rotation  of  the  plate  normal  with  respect  to  the 

crystal  axes  accompanying  a  temperature  change,  which  is  a  result  of  the 

2 

anisotropy  and  was  ignored  in  earlier  work  ,  is  automatically  included 
here.  In  the  description  we  employ  the  changes  in  the  effective  elastic 
and  piezoelectric  constants  have  less  symmetry  than  the  fundamental 
elastic  and  piezoelectric  constants.  As  a  result,  in  the  general  aniso¬ 
tropic  case  there  are  45  independent  changes  in  the  effective  elastic 
constants  and  27  independent  changes  in  the  effective  piezoelectric 
constants  3S  compared  to  21  independent  elastic  constants  and  13  inde¬ 
pendent  piezoelectric  constants.  In  the  case  of  quartz  there  are  ten 
independent  changes  in  the  effective  elastic  constants  and  four  independent 
changes  in  the  effective  piezoelectric  constants,  as  compared  to  six 
independent  elastic  constants  and  two  independent  piezoelectric  constants. 

Since  under  this  program  laterial  field  excitation  is  being  employed 
in  addition  to  thickness  excitation,  no  doubly-rotated  orientations  are 
required  for  the  determination  of  the  temperature  derivatives  of  the 

effective  elastic  constants  even  though  a  larger  number  of  coefficients  is 

2 

to  be  determined  than  heretofore  .  Since  the  piezoelectric  coupling  is 
small  in  quartz,  the  changes  in  the  piezoelectric  3nd  dielectric  constants 
cannot  be  found  with  accuracy  lr„rr.  the  measurement  uf  change  in  thickness 


S 


resonant  frequency  xvith  temperature  of  uncompensated  cuts.  The  changes  in  th 
effective  piezoelectric  and  dielectric  constants  can  be  accurately  deter¬ 


mined  from  the  measurement  of  the  temperature  dependence  of  the  resonant 
frequencies  of  the  fundamental  or  harmonic  overtone  trapped  energy  moces 
in  compensated  cuts,  such  as  AT,  SC  and  BT  cuts.  Clearly,  the  changes  in 
the  ten  effective  elastic  constants  can  readilv  be  determined  from  data  on 
the  temperature  dependence  of  the  thickness  resonant  frequencies  of  the 
uncompensated  cuts.  In  particular,  five  can  be  determined  irom  the 
unrotated  cuts  and  four  can  be  determined  from  the  one ^  rotated  Y-cut 
that  l s  needed  for  the  determination  of  c  ^  Jn<^  c  ^ne  1  cut 

consisting  of  a  rotation  about  2,  preferably  of  about  -o' ,  is  requirou  to 
obtain  the  remaining  one. 

Since  trapping  of  all  modes  will  be  employed  to  eliminate  coupling 

to  unwanted  effects,  all  coefficients  determined  from  the  pure  thickness 

analysis  will  subsequently  be  refined  by  successive  iteration  using  the 

16 

analysis  tor  the  trapped  energy  resonator  .  Finally,  rodunuant  chocxs 
will  be  made  using  other  orientations  and  overtones,  including,  of  course, 
the  thermally  compensated  cuts. 


2 .  Basic  Equations 

The  linear  electroelastic  equations  for  small  fields  suDerposed  on 
a  bias,  which  are  required  in  this  work,  may  be  written  in  the  formw'i"/ 


where 


K,  t  =  2  u  ,  J 
Lv ,  L  o  \  • 


L.L 


=  0 


-  ^ 


n 


f  2.  2) 


q 


and 


K~  =  c  u  4-  p  4  —  p  u  -  c 

Lv  :LvMa  a,M  ’  ~L  LMv  v,M  'LM"\M 


(2.  3) 


\\  AGLvMaU^,M  +  ^Lv^M  ’  ARLMYU-/,M  ' 


(2.4) 


Equations  (2.1)  constitute  the  stress  equations  of  motion  and  charge 

equation  of  electrostatics  referred  to  the  reference  coordinates  of 

material  points  at  the  reference  temperature  T  ,  i.e.,  before  the  static 

o' 

deformation  resulting  from  tne  change  in  temperature  to  T  occurs,  which 


are  called  reference  c 


oordinates  and  are  denoted  .  In  Eqs. 


1 )  and 


(2.2)  Ew  ,  Dt  and  u  denote  the  components  of  the  small  field  Piola- 
Kirchhoff  stress  tensor  which  is  asymmetric,  the  reference  electric 
displacement  vector  and  dynamic  portion  ot  the  mechanical  displacement 
vector,  respectively,  and  denotes  the  reference  mass  density.  In  (2.2) 

'■N.  'V 

for  convenience  we  have  written  both  K,  and  Li  as  the  sum  of  a  linear 

Lv  L 

~x  "n  ~  n 

dynamic  part  and  and  a  nonlinear  static  part  and  dr,  .  The  linear 

dynamic  portions  and  are  the  ordinary  symmetric  mechanical  stress 

tensor  and  electric  displacement  vector  of  linear  Diezoelectricity  and  are 

given  by  the  usual  linear  piezoelectric  constitutive  relations  in  (2,3'*, 

where  L  denotes  the  linear  dynamic  electric  potential.  The  quantities 

ct  w  •  and  e  denote  the  second  order  elastic,  piezoelectric  and 

2  LvMn '  MLv  LM 

dielectric  constants,  respectively.  The  nonlinear  static  quantities 

and  2r^  in  (2.2)  and  (2.4)  are  the  portions  of  the  asymmetric  Piola- 

Kircnhofr  stress  tensor  and  reference  electric  displacement  vector  resultin 

from  the  piezoelectric  vibration  in  the  presence  of  the  biasing  state 

caused  bv  the  temperature  change  (T  -  T  ). 

o 


10 


The  symbols  XG  ,  ^Vn 


and  denote  changes  in  effective 


material  quantities,  which  can  be  expressed  in  terms  of  fundamental 
material  constants  as  shown  in  previous  work^ '  'L '  .  However,  since 
temperature  derivatives  of  the  material  constants  higher  than  the  first 
are  to  be  determined  in  this  work,  it  is  not  feasible  to  try  to  find  the 
temperature  derivatives  of  the  fundamental  material  constants  here  becaut 
the  required  higher  order  fundamental  material  constants  of  quarto  are 
not  presently  known  and  would  be  prohibitively  costly  to  evaluate.  dor.Si 
quently,  we  take  the  alternative  course  of  evaluating  the  changes  in  the 


effective  material  constants  _GT 


~^ML\  •  “  'LM 


in  this  work.  wher< 


V  _L  {T  _  ;n  o~  ■ 

-  L°’  ,_ni 


Before  discussing  the  svmmetries  of  the  GT  w  and  the  XT  and  the 

L\  Met  ML'. 

indicial  notation  employed  in  Eqs.(2.i)  -  <  2.U),  we  consider  it  idvisab  Le 

for  clarity  to  briefly  outline  the  nature  of  the  deformation  that  must  he 

accounted  for  in  the  description.  At  the  reference  temperature  the 

points  of  the  body  are  denoted  by  the  reference  coordinates  X...  '.'Then  the 

temperature  is  changed  from  T  to  T  the  Doints  of  the  bodv  move  to  new 

o 

positions,  which  are  called  intermediate  coordinates  and  are  denoted  by  r  , 
where  f  =  c_  (X_).  Clearly,  the  static  displacement  w  mav  be  denoted 

~~i  ~Zi  L»  ^ 

w  =  ^-X.  When  the  body  is  vibrating  at  some  temperature  T,  the  points  of 
the  body  move  from  to  the  present  position  y.  ,  where  y.  =  y .  •  t '  = 

V, -1  -  'X,  ). t ,  =  y  < X  . t ) .  Clearly,  the  dynamic  displacement  u  mav  be 
denoted  u  =  y  -  \  ana  we  have  y  =  X  *  w  -+■  u.  At  this  point  it  is  also 
purnoseful  to  iote  Cl i.  t  if  at  the  reference  temperature  T  a  noint  n  the 


surface  of  the  body  has  unit  outward  normal  N  relative  to  the  principal 

Li 

axes  of  the  crystal,  when  the  temperature  is  changed  to  T  that  same  point 

has  a  different  unit  outward  normal  v  relative  to  the  principal  axes  of 

or 

the  crystal.  In  the  earlier  work  on  the  temperature  derivatives  of  the 
elastic  constants,  which  was  based  on  the  linear  theory,  the  equations 

could  be  referred  onlv  to  what  are  here  called  the  c  coordinates  and  the 

yy 

difference  between  v  and  N  was  ignored.  Since  all  geometric  measurements 

are  made  at  the  reference  temperature  and  when  the  equations  are 

referred  to  the  coordinates,  the  density,  thickness  ar.d  surface 

normal  N  never  change,  it  is  clearly  significantly  advantageous  to  use 

this  description  Doth  for  simplicity  and  accuracy. 

The  aforementioned  use  of  reference  coordinates  is  the  reason  that 

the  effective  constants  G,  ,,  and  R.,T  have  less  symmetry  than  the  funda- 

L\Ma  MLv 

mental  elastic  and  piezoelectric  constants  cT  w  and  e^,r  .  respectively. 

^  2LvMg-  elLv 

Although  the  cT  admit  interchanges  on  each  pair  of  indices,  the  G.  w 
•>LvMa  l'>  M-, 

have  only  the  symm<  try 


=  S 


mLv  ’ 


(2.5) 


and  have  no  symmetry.  Consequen. ly,  in  the  general  inisotropic  case 

there  are  45  independent  effective  an^  ^7  independent  effective  R^^. 

whereas  there  are  21  independent  ^vMr  and  18  independent  e^.  .  The 
discussion  in  the  above  paragraph  also  makes  clear  the  reason  tor  the 
mixed  notation  of  capital  latin  and  lower  case  greek  indices.  Although 
this  can  be  eliminated  on  the  cr  ..  and  e...  in  Eqs.i2.3),  it  should  not 
be  eliminated  on  the  G^,  and  in  Eqs.(2.4)  because  it  emphasizes 


the  Lack  of  symmetry.  Consequently,  we  retain  the  mixed  notation,  which 


is  consistent  with  the  notation  of  Refs. 4,  14  and  17,  throughout.  The 
cycles  above  variables  have  been  introduced  for  consistency  with  Refs. 14 
and  17.  We  have  employed  Cartesian  tensor  notation  and  vector  notation 
interchangeably  and  the  convention  that  a  comma  followed  by  an  index 
denotes  partial  differentiation  with  respect  to  a  reference  coordinate, 
the  dot  notation  for  differentiation  with  respect  to  time  and  the  summation 
convention  for  repeated  tensor  indices. 

Since  the  quantities  referred  to  as  nonlinear  are  static,  from  (2.1) 
and  (2.2)  we  obtain  the  dynamic  linear  stress  equations  of  motion  and 
charge  equation  of  electrostatics  in  the  form 

"LV.L-  <L-°>  ^7> 

which  with  the  linear  constitutive  relations  (2.3)  yield  the  usual  equations 
of  linear  piezoelectricity.  We  further  need  the  matrices  of  the  elastic, 
piezoelectric  and  dielectric  constants  of  quartz  referred  to  the  princio 

axes  of  the  crystal,  which  is  in  class  32.  These  matrices  may  be  written 

.  -  IS 

in  the  torn 


CLi 

C 1 2 

C13 

°  14 

0 

0 

C 1 2 

'll 

C 13 

"C  14 

0 

0 

C13 

C  13 

°  33 

0 

0 

0 

°14 

'C 14 

0 

C  ,  , 

0 

0 

0 

0 

0 

0 

C44 

C14 

0 

0 

0 

0 

C 14 

°  66 

C66  2  (CU  '  C12)' 


ell 

-ell 

0 

e14 

0 

0 

0 

0 

0 

0 

'e14 

-ell 

0 

0 

0 

0 

0 

0 

J-J 


'll 

0 

0 

0 

c 

'll 

0 

0 

i 

0 

'33 

(2.3) 
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in  which  we  have  employed  Che  conventional  compressed  notation 

The  equation  for  the  perturbation  in  eigenfrequency  obtained  from 

14 

the  perturbation  analysis  may  be  written  in  the  form 


A  =  H  / 2jj  ,  D  =  'll,  -  4,  , 

Ji  ll 


where  u,  and  x>  are  the  unperturbed  and  perturbed  eigenfrequencies , 
respectively,  and 


(2.9) 


H  =  - 


rr,„ 


V 


+  2ARIMvt,Lgy,M  “aLM",  L'.MJ  '  'o 


A'T„fur  r,  dV  . 


(2.10) 


The  vector  g'-  denotes  the  normalized  mechanical  displacement  of  the  j,th 

y 

unperturbed  mode  and  f"  denotes  the  normalized  electric  potential  for  the 
_th  mode,  i . e . , 

~‘/N  , 


u4/n  .  r  =  '?Vn  ,  (2.ii) 


V 


where 


N2  -  J  su^u^dV, 


(2.12) 
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In  (2.11)  and  (2.12)  and  cp^  represent  the  mechanical  displacement  and 

electric  potential,  respectively,  of  the  uth  eigensolution  of  (2.7)  with 

(2.3)  subject  to  the  appropriate  boundary  conditions. 

Since  Eq.(2.10)  contains  ^®MLv  anc*  we  neec*  c^e 

matrices  of  these  quantities  referred  to  the  principal  axes  of  the  quartz. 

Since  c^e  matr^x  f°r  is  c^e  same  as  t^e  matrix  for  given 

in  (2.8),.  However,  before  we  write  the  matrices  for  the  G.  w  and  R  ,  we 
3  mLy  ’ 

must  introduce  a  convention  for  a  larger  range  compressed  notation  than 

the  one  cotnmonlv  emoloved  because  of  the  reduced  svmmetrv  on  G,  and  R... 

l.vMg'  MLv 

compared  with  anc*  resPectavely •  To  this  end  we  introduce  the 

">0 

convention  shown  in  Table  I.  Then,  using  results  of  Mindlin  for  matrices 
having  these  symmetries  for  quartz,  we  obtain  the  required  matrices  in  the 
form 


G11 

G12 

G13 

°14 

0 

0 

G17 

0 

0 

G1 2 

G11 

G13 

-G14 

0 

0 

"G17 

0 

j  i 

i 

G13 

G13 

G33 

0 

0 

0 

0 

0 

0 

G14 

“C14 

0 

G44 

0 

0 

G47 

0 

0 

G 

0 

0 

0 

G__ 

G17 

0 

G47 

G17 

0 

0 

0 

0 

G17 

G66 

0 

G14 

G69 

G17 

-G17 

0 

G47 

0 

0 

G55 

0 

0 

0 

0 

0 

0 

G47 

G14 

0 

G44 

G14 

0 

0 

0 

0 

G17 

G69 

0 

G14 

G66 

L5 


R11 

'Rll 

0 

R14 

0 

0 

R17 

0 

°  ! 

R,,  = 

Mp 

0 

0 

0 

0 

~R17 

-Rll 

0 

-R14 

■Rn  j 

0 

0 

0 

0 

0 

R36 

0 

0 

-R36 

i 

3 .  Pure  Thickness  Vibrations 

We  first  consider  Lateral  excitation  of  thickness  vibrations  because 
in  that  case  the  major  surfaces  of  the  pLate  are  uneiectroded  in  the  /acini ti¬ 
er  the  mode  and,  consequently ,  the  three  waves  are  exactly  uncoupled  at  the 
surfaces  at  resonance  in  the  general  anisotropic  case.  To  this  end  at  the 

reference  temperature  T  we  Locate  the  origin  of  coorainaces  at  the  center 

o 

of  the  pLate  of  thickness  2h ,  with  the  X^-coordinate  directed  along  the 
trigonal  axis  and  the  X?-coordinate  directed  along  a  digonal  axis.  We  denote 

the  unit  normal  to  the  top  major  surface  of  the  Plate  at  T  by  N...  Hie 

o  '  K 

plate  is  driven  into  thickness  vibration  by  an  electric  field  of  magnitude 
i'xt 

E  °  ,  which  at  temperature  V  is  directed  alone  the  unit  vector  S.  in 

o  o  '  x 

the  plane  of  the  plate  and,  hence,  normal  to  M,  so  that  it  satisfies 

K. 


\\  =  0 


(3.1) 


but  is  otherwise  arbitrary.  Now,  substituting  from  (.2.3)  into  (2.7),  we 
obtain 


I  U  T  P  =  ^  u 

LvMa  O', ML  MLv ~ . ML  “o  / 

eLMva-/,ML  '  SLM^.ML  =  °  ’ 


( 3 .  2) 


which  constitute  the  four  coupled  equations  of  linear  piezoelectricitv  in 
u  and  'o  and  where  we  have  taken  the  liberty  of  omitting  the  Lower  script  2 


L  6 


on  c.  .  Since  the  major  surfaces  of  the  plate  are  trac tion- free,  we 
LvMa  ’ 


have  the  boundary  conditions 


"A?  ’ 0 


by 


(3.3) 


and  since  the  electric  field  in  space  vanishes  at  = ±  »  and  we  are 

concerned  with  the  thickness  solution,  we  effectively  have  the  electrical 
boundary  condition 


Vl*  -  ° 


(3.4) 


on  the  major  surfaces  of  the  plate.  The  substitution  of  (2.3)  into  (3.3) 
and  (3.4)  yields 

NL(cLvMorU0f,M  +  eMLv'^,M)  =  °  * 

NL(eMLyUY,M  '  eLMCp,M)  =  °  ’  3C  W  =  ±h  *  (3-5) 

As  a  solution  satisfying  the  differential  equations  (3.2)  and 

o  L'-ut 

boundary  conditions  (3.5)  with  a  driving  electric  field  E  S„e  ,  we  take 


«,-*v  -  W1"'  5*!B  sin'VlC+E“<SKXK‘r,W!  »•*> 


10Jt 


where 


r  “  NK£KMSM//NReRSNS  ’ 


(3.7) 


and  is  required  to  assure  satisfaction  of  the  boundary  condition  (3.5)„. 
The  solution  in  (3.6)  satisfies  (3.2)  provided 


hlV 


c  5  )  A  -  0  , 

yet  ot  ’ 

(3.3) 

:lwyAy/nknr£kr  » 

(3.9) 

where 


L  7 


(3.10) 


c 

ya 


ya 


+  e  e  /e 
y  at 


2  2 

c  =  o 

o 


and  we  have  Introduced  the  notation 


Cva  WlvMQ  ’  ey  NRNS6RSy  »  S  ^L£LM^M  ’  (3.11) 


which  is  convenient  in  this  work  and  has  nothing  to  do  with  the  compressed 
(p,q)  notation  in  Section  2.  Equations  (3.S)  constitute  three  homogeneous 
linear  algebraic  equations  in  the  and  the  condition  for  a  nontrivial 
solution  yields 


c  -  cS  =  0  . 

vor  var 1 


(3.12) 


Equation  (3.12)  is  a  cubic  in  Z,  which  for  a  given  yields  three  real 
positive-’*’  c(n)  (n  =  1  2,  3),  which  we  assume  to  be  distinct.  Each 
yields  amplitude  ratios 


(n)  (n)  (n) 

1  :  A2  :  A3  ’ 


(3. 13) 


when  substituted  in  any  two  of  the  three  equations  in  (3.8),  where  the  A_ 
satisfy  the  orthogonality  relations 


A(’*)A(V> 


(3.14) 


and  N ^  is 
ratios  by 


the  normalization  factor. 


Defining  the  normalized  amplitude 


Q(^> 

y 


<3.153 


we  obtain 
functions 


the  orthogonal  matrix 
(3. 6)  with  (3. 7)  into 


5  .  The  substitution  of  the 

jot 

the  boundary  conditions  (3.5) 


solution 

yields 
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-  0  . 


(3.16) 


~1  c  A  cos  "  h  +  a/E° 
ya  y  y 

and  (3.5),,  is  satisfied  identically  by  virtue  of  (3.9),  where 


/ 

e 

v 


NLSM6MLv 


rIW 


MLv  ‘ 


(3.17) 


From  (3.16)  we  see  that  the  exact  condition  for  pure  thickness  resonance 
under  lateral  field  excitation  is 


cos  '  h  =  0  or  a  =  nr/  2,  a3  1,3,5  .... 


;  5.18' 


With  the  aid  of  the  orthogonality  of  the  5  we  mav  write  ■'  3.  lb)  in 

a  particularly  illuminating  form,  in  which  each  thickness  noae  is  uncoupled, 

simply  by  referring  the  equations  to  the  coordinate  system  consisting  of 

the  eigenvector  triad.  Although  the  form  is  not  of  great  use  to  us  in 

this  case  of  lateral  excitation,  it  is  important  in  the  case  of  thickness 

excitation,  which  is  treated  next.  To  this  end  we  transform  (3.16)  to 

the  eigenvector  triad  with  3  and  write 

w 


A  -  3  A  A  =  3A 

mot  or  or  |ua  q,  7 

and  substitute  from  (3.19),,  into  (3.16)  and  employ 
eigenvectors  3^^  and  make  use  of  the  orthogonality 


(3.19) 

(3.3)  for  the  normalized 
of  the  3  to  obtain 


cos 


h  +  3 


\a  v 


=  0  . 


:o> 


which  constitute  three  uncoupled  equations  (one  for  each  v)  giving  the 

~  o 

amplitude  A  in  terms  of  the  driving  amplitude  E  . 

In  the  case  of  thickness  excitation,  which  we  now  treat,  we  still 
have  the  differential  equations  (3.2)  and  the  boundary  conditions  (3.31. 
'out  instead  of  ' '  3  .  V) .  we  have 

19 


y  =  i  7  e^*'  at  =  ±h  .  (3.21; 

On  account  of  this,  although  the  solution  for  u^  is  still  of  the  form  shown 
in  (3.6)  the  solution  for  'a  now  takes  the  form 

?’[B  sin  -  \\+(c  +|)  Y\]  «l”t-  '  — > 

These  solution  functions  satisfy  the  same  differencial  equations,  i.a., 
(3.2),  as  in  the  case  of  lateral  excitation  because  the  expression  for  s 
in  (3.22)  differs  from  the  one  in  (3.6),  only  by  terms  linear  in  >2, . 
Consequently,  £qs.(3.8)  -  (3.15)  still  hold  for  the  case  of  thickness 

excitation.  Now,  substituting  from  >3. 22)  into  ;'2.21)  ana  emo loving  13.  Li) 
we  obtain 

C=-(eA/eh)sin~h.  (3.  23) 

y  y 

As  noted  in  the  Introduction,  for  the  case  of  thickness  excitation  the 
analysis  is  restricted  to  orientations  for  which  the  three  coupled  waves 
at  the  conducting  surfaces  of  the  plate  essentially  uncouDle  and  one 
approximately,  but  very  accurately,  dominates  each  uncoupled  vibration. 
Under  these  circumstances  only  that  component  of  the  traction  vector  in 
the  direction  of  the  eigenvector  of  the  dominant  wave  need  be  considered 
at  a  time.  The  successive  cons ideratiou  of  the  three  distinct  eigen- 
directions  yields  the  three  approximate,  but  very  accurate,  uncoupled 
solutions.  To  this  end  we  substitute  from  (3.6),,  (3.22)  and  ‘3.23) 
into  (3.5).,  which  is  equivalent  to  <3. 3),  and  transform  13.5).  with  5, 

L  L  VV 

and  substitute  from  (3.19),  and  employ  (3.8)  for  the  normalized  eigen- 

(  ) 

vectors  6  ~  and  the  orthogonality  of  the  5  to  obtain 
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where 


cos 


A  /V 

e  e 
'J  v 


sin 


'N, 

eh 


v  2 


(3.24) 


e  =  S  e  . 

v  w  y 


(3. 25) 


Equations  (3.24)  are  the  uncoupled  equations,  one  for  each  v  <v=  1,2,3), 
which  give  the  amplitude  A  for  the  vth  pure  thickness  mode  in  terms  of  the 
driving  voltage  V.  From  (3.24)  we  obtain  the  condition  for  resonance  of 
the  vth  mode  in  the  form 

Can  (v)h  =  (v)h/kv  ’  ' 3. 26 1 

where 

k~  =  (ev)2/c(v)?'.  (3.27) 

The  forms  shown  in  Eqs.(3.24)  and  (3.26)  were  obtained  because  for 

convenience  the  mass  loading  due  to  the  finite  thickness  of  the  electrodes 

was  ignored  in  Eqs.(3.3)  and  13.5)^.  If  the  mass  loading  had  been  included. 

22  °3 

in  place  of  (3.26)  we  would  have  obtained- 


tan  ,h  = 
(v) 


■,  .h/(k2+RT“  . h ~)  . 
(v)  V  (v) 


13. 23) 


where  R  =  2z'h' / z  h  and  z'  and  2h'denote  the  mass  density  and  thickness  o! 
coo 

an  electrode,  respectively.  Since  the  piezoelectric  coupling  is  small  in 
quartz,  the  usual  expansion  of  (3.29)  about  "  ^h  =  nn/2  (n  =  1.3.5,  .  .  .  ^ 


with  the  aid  of  (3.10)-,  yields' 


•u 

n 


,-(v)x  1/2. 

=  sn  (s _ ^ 

2h  \  o  ,  v 


n 

4k“ 


_ v_ 

?  -> 


(3.30) 
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from  which  we  readily  obtain 


O-  -  3w  )/3oi  =  32  k~/9rr  , 

3  1  o  v 

2 

to  order  k“  and  where  for  small  coupling 


(3.31) 


0) 

o 


(n/2h)(c(v)/o  )1/2 
o 


c<v)  =  c  <v)d-k2/2)  . 

V 


(3.32) 


Equation  (3.31),  with  (3.32)  is  useful  for  obtaining  the  piezoelectric 
constants  from  measurements  of  the  fundamental  and  third  overtone  thicknes 
excited  thickness  resonances. 

In  order  to  use  the  foregoing  general  anisotropic  results  for  pure 

thickness  resonances  to  evaluate  the  material  constants  from  data  at 

T  =  25“C,  we  must  consider  various  specific  orientations  of  quartz .  To 
o  ’ 

24 

this  end  we  introduce  the  conventional  IEEE  notation  for  doubly-rotated 
cuts  of  quartz  and  write  (Y,  X,  w,  i)co,  9,  which  yields  the  relations 


s 


M.  =  -  cos  5  sin  c,  N_  =  cos  9  cos  N„  =  sin  9  ,  Cl.  33) 

1  2  _> 

for  the  components  of  the  normal  to  the  major  surfaces  of  the  plate  at 

in  terms  of  the  rotation  angles  cs  and  9.  In  particular,  we  need  the 

expression  for  c  .  e  .  s.  e  and  c  for  doublv-rota ted  cuts  of  quartz. 

vi’  y  y  YO' 

With  the  aid  of  (3.11),  (3.17)  and  (3.10)^  and  the  matrices  in  (2. S'),  we 
obtain 

C11  =  N1N1C11  +  N2N2C66  +  2N2N3C14  +  N3N3C44  •’ 

=  12  =  N1N2(C12+C66}  +  2N1N3C14-  S3  =  +  WC13 +=44>  ’ 

C 22  =  N1N1C66  +  N2N2C11  "  2N2N3=14  +  N3N3=44  ’ 

=  23  3  Vl=14  *  N2N2C14+N2N3<=13+C44)* 

c.  =  N  N.c  .  +  N  N  c  ,  -t-N,N0c  .  (3.34) 

3  3  L1-4.U  11  4h  j33-> 


1  *1 


The  substitution  of  (3.43)  into  (3.8)  yields  the  two  uncoupled  systems 


!c,,  -c)  rv,  ~0,  i.cA^  -  c  '  A0  4  c  -  3 .  -liLA,  tc 


dd 


'14* 


-  a 

44  "3 


the  first  or  which  vie  Ids 


_■  1) 


'11  '11  1 


with  amplitudes 


a!"  ’  =  'I,  . 

a 


).-o> 


Touaticns  ■  3 . — *  '<  ..  ,  yielu 

^dd  -i.4  — i  nd  14 

which  is  a  quadratic  equation  yielding  the  two  roots  o'  and  c  '  with 
amplitude  ratios 


From  Eg.  3.4  7')  we  obtain  the  two  relations 


_! 2)  ^  _(3) 

C  -r  C 


-( 2)_<  3) 


i  3.  ^9' 


which  prove  to  hie  use  tul  in  tnis  work.  liquations  •  1.  +2)  reveal  that  :  >r 
the  X-cut  plate  we  have  thickness  (X)  excitation  or  the  p lecoe  loc tr ica i 1 v 
stiffened  extensional  mode  and  lateral  'V'  excitation  u  both  purely  ’  iasti. 
shear  modes  which  are  couol«?u. 
c o r  ~he  Y-cat  we  have 


XI,,  U,.  N^)  =  <0,1.03  . 


which  with  (3.34)  -  (3.37)  and  (3.39)  yields 


-12  C 1 3  ~  C22  Cll-  °23  "  °  14’  °33 


'1  11’  2 


=  0. 


11 


’  el  °’  e 2  *  Siell’  e3 


1  14 


’ll  =  °66  ""311  ’311’  ?12  °’  °  13  322  CU’  r:3  C14’  C33  Cw4 


The  substitution  of  (3.53)  into  (3.8)  yields  the  cvo  uncoupled  ;vst. 


:  )  A1  0  .  fcu  -  c)  A  ~  0  • 


-  ci^A  ,  "  (ci,  '  A3  "  0 


the  first  of  which  yields 


_  (  l )  ^ 

C66  °11  "'ll 


with  amplitudes 


(1.0 


i 


Eauations  (3.54),  ,  yield 
_ .  j  ' 


c  -  {c 


11 


=  0 


_ . : ) 

which  is  a  quadratic  .aquation  yielding  the  two  riots  c  mu  c 
amplitude  ratios 


i  n  i 


1  A, 


1  n  i 

)  a 


x,n)i  = 
3 


). 


.  i  n  • . 


14 


"rom  Ea . (  3.57)  we  obtain  the  two  relations 


2)_v 3) 


11 


11' 


which  are  useful  in  chis  work.  Equations  (3.52)  reveal  than  for  the  Y-cut 
plate  we  have  thickness  (Y)  excitation  of  the  piezoelectrical lv  stiffened 
(X)  shear  mode  and  lateral  (X)  excitation  of  the  other  two  purely  elastic 
coupled  shear  and  extensional  modes. 

For  the  Z-cut  we  have 


(Np  M3)  =  (0,0.1)  .  O.  60) 

which  with  (3.34)  -  i3.37)  and  (3.39')  yields 


C1L  Ci4-  "L2_0-  c13  °-  c;2°c44'  c23'°'  "33  c33 


o.  s  -  o,  a  -  4  j,  .{-0.  *3  *  ° 


C  -  C 


<  3. hi) 


(  3. 62) 

(3.63) 


vcr  V  ot 

The  substitution  of  (3.63)  with  (3.61)  into  (3.8)  yields  the  three  uncoupled 
equations 


(c , 


Al  =  °. 


(  c 


44 


-  c)  A t  =  0 > 


•s = ° 


which  yield 


-d) 


-(3) 

c 


=  c 


33 


along  vith  the  matrix  of  amplitude  ratios 


(  3  .  b4 ) 


(3.65) 


(3 , ob) 
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Equations  (3.62)  reveal  that  for  the  2-cut  plate  we  do  not  have  thickness 


excitation  of  any  modes  and  we  have  lateral  (X)  excitation  of  Y-shear  onlv. 

Since  the  analysis  for  the  rotated  X-cut  yields  no  significant 
simplification  over  the  analysis  for  the  doubly-rotated  plate,  we  do  not 
present  che  results  for  this  case  here.  We  simply  note  that  all  three 
coupled  modes  can  be  driven  by  both  thickness  and  lateral  excitations. 

For  the  rotated  Y-eut  we  have 


N,  ,  N\.  N, )  =  (0.  cos  9,  sin  9) 
»’hich  with  (3.34)  -  (3 .37)  ana  i'3. 39>  yields 

=  ■) .  i ,  .  =  0 . 


J  .  0/ 


■* 

■) 

- 

c 

11  =  C  “66 

-  -CSC.  ,  -1-  3  “c 

i4‘  “12 

•> 

“> 

c 

n  "  c  '-oo 

-  2csc,  .  +  s “c 
14 

->■-+  -3 

0 

C  33 

C "C44  * 

- 

n 

e ,  =  -  o 

>ii  1  '  ^5  C  ,  ,  , 

e  ,  =  0 .  - 

- 

:I  l  *4 

- 

/ 

.  f 

,0 

wl 

=  ‘  “Y°!i 

•> 

=  “  Ve: 

=  -4-  - 

"13  =  0- 

'll 

"11 

~  ’  ~12  =  0; 

O 

which  c 

and  s  raDresent  cos  8 

and  sin 

13 


s  “c 


33' 


11  '1 


;e .  .  .  a. 


3.  68) 


s  -,0 


(3.70) 


of  '3.70)  into  (3.3)  yields  the  two  uncoupled  systems 


uii'A)  Aisl)'  A:  ~C:3A3  “ 0 

3  51  V’ 


the 


;t  ot  which  vie 


las 


_  !  i  ) 


3.  72) 


with  amplitudes 


which  are  useful  in  this  work.  Equations  (3.69)  reveal  that  for  the 
rotated  Y-cut  plate  we  have  thickness  (rotate 3  Y)  excitation  of  the 
piezoelectricaily  stiffened  <, X )  thickness-shear  mode  ana  lateral 
>.!0  excitation  of  the  other  two  purely  elastic  coupled  shear  anu 
extensional  modes. 

From  the  foregoing  it  is  clear  that  all  the  elastic  constants 

can  be  determined  from  measurements  using  lateral  excitation.  First 

c,.  can  be  determined  from  a  3-cut,  then  c,  ,  can  be  determined  from 
-*4  11 

a  Y-cut  and  c..  from  an  X-cut  3nd  c . .  from  either  a  Y-cut  or  an  X-cut. 
ob  14 

Finally,  both  c,.,  and  c„.,  can  be  determined  from  one  rotated  Y-cut. 

It  is  also  clear  that  the  piezoelectric  and  iielectric  constants  can 
be  determined  from  thickness  excitation  of  the  fundamental  and  third 
harmonic  of  a  Y-cut  or  m  X-cut  and  three  rotated  Y-cuts. 


We  now  present  some  of  the  results  for  a  rotation  about  Z 

even  though  they  are  not  needed  for  the  determination  of  the  constants 

at  T  =  25°C  because  one  is  needed  for  the  determination  of  the 
o 

changes  in  one  of  the  effective  elastic  constants.  For  a  rotation 
about  Z  we  have 


(Np  N.,,  N3)  =  (-sin  cp,  cos^,0)  , 
which  with  (3.34')-  (3.37)  and  (3.39)  yields 


(3.77) 


C11  ~  S  C11+'C  C  66  ’  C12  ‘SC(C12+C66)’ 


c  13 


23 


2sc  ,  c 

22 

S  C66  _L  C  C11  » 

2 

0  o 

14  *  C  C 14- 

C33 

=  s~c,  .  +  c  "c 

*44  4-4 

=  c 


2  2-.  -  2  2 
s  e  i -  c  e,  .  .  e,  =  2sc  e. . ,  a  =  s  a,,  +c  e 


(3.73) 


11  '  "  C11  "  "11  ’ 


e]_  “  ”  ^oc  ”  ^r,cel4  ’ 


e 2  ~  ‘  S1C  ell  +  S2Sell’  e3  =  S1C  el4  *  b2Se14  ’ 


(3. 79) 


‘  1 


_  -  ^  .  1  2  -  _  ^ 

f  —  r*  f  -  c  =  c  -*■  ■—  r  =  r 

11  11  s  ’  c12  12  e  L1  »  13  13’ 


=  c 

00  00 
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C  =  C  C  =  c 

23  23  ’  33  33  * 


(3.30) 


Since  the  substitution  of  (3.80)  into  (3.8)  yields  no  significant 
simplification  over  the  doubly-rotated  plate,  we  do  not  present 
the  detailed  expressions  for  the  coupled  linear  equations,  the  3  •  3 
determinant  or  the  amplitude  ratios  for  the  plate  rotated  about  Z. 
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4 .  Perturbation  of  Thickness  Vibrations 

We  first  obtain  the  one-dimensional  expression  for  the 

perturbation  integral  H  from  the  general  three-dimensional  expression 

given  in  Eq.(2.10)  because  the  one-dimensional  version  is  directly 

useful  for  the  treatment  of  thickness  vibrations  discussed  in 

Section  3.  To  this  end  we  write  dV  =  A  ds  in  Eqs.(2.10)  and  (2.12) 

o  o  o 

and  then,  since  the  pure  thickness  solution  does  not  vary  along  the 
surface  of  the  plate,  factor  out  the  A  from  both  expressions,  which 
cancels  out  of  the  entire  description  by  virtue  of  (2.11).  Under 
these  circumstances  the  expressions  for  H  and  N  take  the  respective 
forms 


H  = 


IT- 

v  L 


-h 


GLvMoSa,MSv,  L  +  2ARim 


M 


f“.  !  ds 


=LM  ,  L  ,MJ 


4.1) 


-> 

N“  = 


r  ds 


-h 


v  v 


(4.2) 


where 


s 

o 


=  Vk  • 


(4.3) 


Bv  virtue  of  the  one-dimensional  [thickness  (s  )]  dependence,  with 

o 

the  aid  of  the  chain  rule  of  differentiation,  we  may  write 


u  M  =  (au  /ds  )(cs  /*x  )  =  N  5u  /5s  .  14.4) 

a.M  arooMM^o 

From  (3.6),  (3.22)  and  (4.3)  we  see  that  in  the  present  notation  the 
thickness  eigensolutions .  i.e.,  with  V  =  0.  can  be  written  in  the  form 


u  =  A  sin  Is  ,  j=B  sin  's  +  Cs  , 
v  v  o’  o  o  ’ 


(4.5) 
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where  B  is  given  by  (3.9)  and  C  is  given  by  (3.23)  for  Che  case  of 

thickness  excitation  and  C  =  0  for  the  case  of  lateral  excitation 

,  .  Lxt 

and  we  nave  suppressed  tne  e 

Substituting  from  (4.5)  into  ('4.2)  and  performing  the  inte¬ 
gration.  we  obtain 

N~  -  c  hA^A"  .  (4.6) 

u  o  y  v 


Now,  let 


A  =  a'-/N  ,  3^  =  3^/N  .  C"  =  C"/N 


then  rrom  (2.11),  (4.4),  (4.5)  and  (4.7),  we  obtain 


,M  =  V'  cos  "So  ’ 


\  =  Nt  (B^~  cos  "s 

L  L  ? 


Now,  substituting  from  (4.8)  into  (4.1),  performing  the  integrations 
and  employing  (4.7),  (3.9),  (3.11)?  and  (3.23),  we  obtain 


a"  -  -  Ac  -Ji»  4  -  o  a'a-A  -  {.& 

-  1_  \  3L  V  J  V  a  „ 


-1 

~y  "\‘ ' 


h 

k9) 


for  the  electroded  thickness  excited  plate,  where 

and  since  C"  =  0  for  the  unelectroded  laterally  excited  plate,  the 
second  expression  in  brackets  vanishes,  we  obtain 


10) 


H 


'J 


(4.11) 


for  the  unelectroded  laterally  excited  plate. 

In  order  to  use  the  foregoing  general  anisotropic  results  for 

thickness  resonators  to  evaluate  the  temperature  derivatives  of  the 

effective  material  constants  CT  ..  ,  and  '  „  from  data,  we  must 

LyMa'  1MY  •'LM 
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consider  various  specific  orientations  of  quartz  as  in  Section  3. 


In  particular.  :;e  need  the  expressions  for  AG  ,  AR  and  A I  for 

'.'OC1  V 

doubly-rotated  cuts  of  quartz.  With  the  aid  of  (4.10)  and  the 
matrices  in  (2.13)  and  (2.8K,  we  obtain 


A, 

4Gn 

=  X1N1AGU 

+  N2N2AG66  +  2N2N3AG17  +  N3N3AG55  ’ 

J.  — 

=  M  N  AG 

12  12 

+  "lVG69  +  2N1N3AG17  • 

1G12 

=  2N  N  AG. . 

12  1-4 

+  WG13  X  N1N3-G2i7  > 

—^2  2 

=  N1N1“G66 

+  N2N2AGil  *  :N2N3AG17  ’  'l3S3AG55  * 

-v. 

“G23 

=  VlAG14 

-  m2n,ag14  -  m2n3ag13  -  s2n3ag.7  . 

“G33 

=  N1N1AG. . 

+  N2N2AG44  +  Wc33  • 

“G66  +  “G69 

=  AGll  "  AG12  ’ 

(4. 12) 

2Rl 

=  n1n1aru  - 

n2n:ar11  -  n,n3ar17  -  n.x_r36. 

AR„ 

=  -  ''N  N  wR 
i  2  11 

*  Wh;  *  •  -R1 ' 0  • 

4.  13) 

a;  =  N,  NT 

i  1 

A  r  +  M  M  a  r  -4-  m  m  A  r 

~bil  ■’  .n  1  3  3  ~33  ■ 

'4.14) 

which  with  the  results  of  Section  3  readily  enables  the  determination 
of  the  relations  we  need  for  any  orientation  we  wish. 

For  X-cut  quartz  we  have  (3.40),  which  with  (4.12)-  (4.14) 

yields 


2.G,  , 

-r-4 


(4. 13) 


(4. lb) 


The  substitution  of  (4.15)  and  (4.16)  along  with  (3.46)  into  (4.9) 
for  the  thickness  (X)  excited  extensional  mode  yields 


'll 


4GU  +  2ARU 


■¥}  W 

c* 

“11 

-  f2iRn  ITT  +  iC 


2  -1-1 
A,  A, 


Ill  +  Ar  11]  11 
11  m 


(4.17) 


The  substitution  of  (4.15)  and  (4.16)  along  with  (3.48)  into  (4.11) 
for  each  of  the  laterally  excited  two  coupled  shear  modes  yields 


u  r 
« t  -I 


A *A^  +  AG, ,  A’A  =1  '''I  h  ,  t  =  2 ,  3  . 


“G66A2A2  +  2AG14“2"3  '  U',44"3“3J 


(4.13) 


For  Y-cut  quartz  we  have  (3.50),  which  with  (4,12)-  ('4.14) 


yields 


-G11  =  AG66-  -G12  =  0-  “G13  ~  ° ’  “*G22  =  "Gll’  ~G23  =  *  AG14'  AG33  =  AG44  • 

(4. 19) 

-Rl  =  '  4RU.  4R:  =0.  4C  =  4CU  .  (4.20) 

The  substitution  of  (4.19)  and  (4.20)  along  with  (3.56)  into  (4.9) 
for  the  thickness  (Y )  excited  (X)  shear  mode  yields 

t  t  ■”  _  '  l  o  a  r»  -  _  \  -  _ 

“1  '  L_“°66  '  ““"11  e,  ,  “*11 

'"'ll 

2  -1-1 
eul  A1 

(4. 21) 


11  ell~  t-1oU2.  T  fu 

ii  11  111  *-  11  sn 


_ni  11 

2  J  h 
ell 


The  substitution  of  (4.19)  and  (4.20)  along  with  (3.58)  into  (4.11) 
for  each  of  the  laterally  excited  other  two  coupled  modes  yields 


=  -  [-11  41 ' 2AG14A11  +  AG44A3A3 


"  5-2,. 


(4.22) 
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For  Z-cut  quartz  we  have  (3.60'),  which  with  (4,12)  -  (4.14) 

yields 


AG 


55’ 


=  AG 


55’ 


AG,^  =  0,  AG^^  =  AG , 

(4.23) 


a  r 


ac33  . 


(4. 24) 


Since  for  the  2-cut  we  have  (X)  excitation  of  (Y)  shear  only,  the 
substitution  of  (4.23)  and  (4.24)  along  with  (3.66)  into  (4.11) 
yields 


H 


U 


= -  ag55a2a2 


(4.25) 


As  in  Section  3  and  for  the  same  reason,  we  do  not  specialize 
the  results  fof  the  rotated  X-cut.  For  the  rotated  Y-cut  we  have 
(3.67),  which  with  (4.12)  -  (4.14)  yields 


AGn  =  c2AG^  +  2scAG17+s2AG,„  AG,,«0,  AG,  ,  =  0  , 


J13 


“Gll  '  2csA^17  AG55, 

AG23  =  -  c'AG14 

■> 

2 

+  csaG^y.  AG22  =  c  AG^ 

-rs“AG33  , 

(4.26) 

7 

~ARU  -  csAR17,  AR2  =  0, 

a;  = c  a:u  +s  a;33 . 

(4.27) 

The  substitution  of  (4.26)  and  (4.27)  along  with  (3.73)  into  (4.9) 
for  the  thickness  (rotated  Y)  excitation  of  the  (X)  shear  mode  yields 

-2 

fi*,  ■ 

■*’=  -A 


H 


I  --[-su  +  2ih 


-2  -1-1 
e,  A ,  A , 


1  11  -1-1  ’  r  -  l  -  cir 

—  -  AC  77J  aiai  P1  -  2[*2AR,  —  +  a:  -4  1 


(4. 231 

The  substitution  of  (4.26)  and  (4.2/)  along  with  (3.75)  into  (4.11) 
for  each  of  the  laterally  excited  other  two  coupled  modes  yields 
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[J  (  --v  '■vC  '•v.  A.r/vr 

H,  -  -  !  +  2“G^3A2A3  +  AG33A|A3 j  "7n  ,  '=2.3. 


<4. 29) 


yields 


For  a  rotation  about  Z  we  have  (3.77),  which  with  (4.12)  -  (4.14) 


^11  =  S'-Gll  C ~AG66’  “G1 2  =  '  SCAG12  *  SCAG69  ’ 

^  ^  ?  o 

— G., =  ~  -Sc4G.  ,  ,  ~  s  1G, .  +  C  2.G.  ,  . 


— Gn„  -  s  4G,  -  c  2.G.  ,  ,  _G„ „  =  s  J.G,  ,  ~c  2.G.  ,  -  _G.  ,  , 

2o  la  14’  jj  j4  -+4  — i 


AR-  =  s  ARu-c_3Ri1.  4R_  =  2sc4R.  , ,  42=  . 


4.  JO) 


,31) 


Since  the  substitution  of  (4.30)  and  (4.31)  into  '’4. 9)  a: 


■ .  _  1  i 


yields  no  significant  simplification  over  the  doub ly-rota tec  place, 

e  U 

we  do  not  present  the  detailed  expressions  for  H  and  H  . 

From  the  foregoing  it  is  clear  that  all  10  independent  _G. 
can  be  determined  from  measurements  of  the  temperature  dependence  of 
the  resonant  frequencies  of  thickness  vibration  of  the  orientations 
considered  in  this  section.  In  particular  £G,  .  .  4G,  ,  .  2.G.  ,  iG,. 


11’  ^  1( 


-a  00 


and  4G._  can  be  determined  from  measurements  of  the  three  anrotated 
oo 

cuts.  In  addition,  4Gn_,  4G,  2.G.  _  and  4G,-  can  be  determined  from 

1  /  a  /  13  jj 

the  one  rotated  Y-cut2^  that  is  required  for  the  determination  of 

c,,  and  c ^ .  Then  4G, „  can  be  determined  from  one  cut  consisting  of 
LJ  jJ  12 

a  rotation  about  Z  and  finally  ±GCC.  can  be  obtained  from  <4.12)-., 

by  / 

which  can  also  be  used  as  a  consistency  condition  for  the  measurements. 


Furthermore,  AR^  and  £R  can,  in  principle,  he  obtained, 

along  with  and  ,  from  measurements  of  the  temperature 

iependence  of  the  resonant  frequencies  of  the  thickness  vibrations 
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of  specific  quartz  plates,  including  at  least  one  doubly-rotated 


cut.  which  is  required  for  However,  in  practice  this  is  not 

possible  on  account  of  the  small  piezoelectric  coupling  in  quartz, 

which  for  most  cuts  causes  the  influence  of  these  coefficients  to  be 

masked  bv  the  influence  of  the  AGT  w  .  Nevertheless,  these  quanti- 

u'/Ma 

ties  can  be  obtained  from  measurements  of  thermally  compensated  cuts 

such  as  the  AT,  BT  and  SC  cuts  (and  others)  because  for  these  cuts 

the  influence  of  the  oertinent  AGT  vanishes.  This  is  discussed 

LvMa 

in  the  next  section  because  the  analysis  including  trapping  is 
required  on  account  of  the  accuracy  needed  to  calculate  these  small 
auantities.  Einallv.  ic  should  be  noted  from  (d.13)  that  _R,  ,  cannot 
be  obtained  from  measurements  of  thickness  resonance  or  have  any 
influence  on  the  behavior  of  any  resonator  vibrating  in  an  essentially 
thickness  mode. 

5 .  Trapped  Enerny  Resonator 

In  this  section  we  briefly  discuss  the  solution  for  the  tr.o.^ec 

.  .  .  .  , .  -  .  ,  1  o .  7 

energy  resonator  witn  significant  reterence  to  previous  work  ' 

because  of  its  length.  However,  before  we  do  this  we  must  aiscuss 

che  fact  that  the  solution  for  the  trapped  energy  resonator  is 

referred  to  the  plate  axes  as  well  as  the  orthogonal  axes  of  the 

eigenvector  triad  of  the  pure  thickness  soLution,  while  in  the 

previous  sections  of  this  work  the  coefficients  are  referred  tu  the 

principal  ixes  of  the  cryscal.  When  the  conventional  IEEE  notation 

discussed  earlier  is  employed,  the  rotation  angles  't  ana  3  ire  the 

first  two  Euler  angles  which  determine  the  orthogonal  transformation  a 

wU 

from  the  crystal  axes  to  the  plate  axes,  where  the  a_  are  given  bv 
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a 


3G 


COS  " 

-cos  &  sin  <p 
sin  3  sin  cp 


sin  c  0 

cos  G  cos  cp  sin  6 

-si;.  3  cos  9  cos  S 


5.  n 


Then  Che  transformation  relations  for  all  the  material  tensors 
employed  in  this  work  are  given  by 


eKNS 


2  K.pN6  ~  aKLar/aMMa6aCLvMa  •’ 

/  __ 

aKLaNMa3veLMv  ’  "'KN  ~  aKL3NM€LM  ’ 

G '  =  a  a  a  a  G 

K3N6  KL  3'/  NM  5a  LvMa 

aKLaNMa3vRLMv  *  CKN  “  3KL3NM3LM  •' 


o.i'i 


in  which  the  primed  quantities  are  referred  to  the  plate  axes.  In 

addition,  since  the  displacement  field  is  referred  to  the  orthogonal 

coordinates  of  the  eigenvector  triad,  which  is  given  by  the  orthogonal 

transformation  3  between  the  orincioal  •’xcs  of  the  crystal  and  the 
VOf 

eigenvector  triad  in  the  preceding  sections  of  this  work  ana  is  given 
by  the  orthogonal  transformation  Q  between  the  plate  axes  and  the 
eigenvector  triad  in  Ref. 16,  we  have 

3  =  Q  a  '5.3* 

\a  w  yc* 

In  this  section  we  use  the  coordinate  convention  of  Ref. 16,  whereby  the 

s  axis  defined  earlier  in  this  work  is  in  the  direction  of  the  X-- 
o 

coordinate  axis  of  this  section  and  X^  and  X^  lie  in  the  plane  of  the 
plate.  We  are  now  in  a  position  to  discuss  the  solution  for  the 
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trapped  energy  resonator  in  relation  to  the  measurement  of  the 
cons  tants . 

It  is  shown  in  Ref. 16  that  the  homogeneous  form  of  the  differ¬ 
ential  equation  governing  the  transverse  behavior  of  the  nth  odd 
harmonic  familv  of  modes  is  given  by 


M 


‘1 


sx: 


n  *1*3 


1  „  U1  n“~“  ~(l)~n  ”n 

p  - _  .  - —  c  u  =  :u 


4h' 


1  ' 


(5.4> 


where 


^  =  u^1  «Xr  X  t)  sin  (nrX,/2h)  . 


o .  ;  i 


and  M  ,  Q  and  P  are  given  in  Eqs.(74)  of  Ref. 16,  and  from  (78),  of 
n-  n  n 

Ref.  16,  we  have 


c<l>  =c(1)(l-8k“/n-2-2R). 


(5.6) 


and 


Iz'h'/z  h  , 
o  o 


<5.7  i 


in  the  notation  of  this  work,  where  z'  and  2h '  denote  the  mass  densit 

o 

and  thickness  of  an  electrode.  Equation  (5.4)  is  for  the  electroded 

region  of  the  plate  and  holds  for  the  unelectroded  region  provided 

- (1)  .  ,  .  .  _( 1)  .  . 

c  is  replaced  by  c  .  If,  as  ls  sometimes  the  case  m  this  worx, 

the  trapping  is  produced  by  an  insulating  film  instead  of  an  electrod 

8(1)  in  (5.4)  is  replaced  by  c*^  where 


:  U) 

C 


_(i) 


(1  -  2R)  . 


(5.8) 


and  R  is  for  the  insulating  film  rather  than  an  electrode.  We  now 

take  the  arbitrary  course  of  neglecting  Q  in  (5.4)  in  order  to 

n 
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eliminate  the  mixed  derivative  term 


Tins 


:n  ora..-r 


is  necessarv 

permit  a  convenient  reoresenta tion  of  the  trapped  merry  mod.1  :  or  an 
arbitrary  orientation  of  the  rectangular  electrodes  or  insu  Latina 
films,  which  are  being  employed  in  the  experiments  beinr  performed. 
It  can  be  shown  that  this  procedure  does  not  result  in  lDpricraoie 
error  except  in  certain  unusual  cases  that  will  not  arise  here. 

It  has  been  shown  in  Sect. Ill  of  hef.lt  that  for  '  -tanruiar 
electrodes  of  Length  1.1  along  X.  and  lb  along  X„ ,  tne  tunuamenta. 
trapped  energy  eigen-solution  for  iav  n  :an  be  written  *r.  the  arm 


where  ,  S,  T  and  C  represent  the  eloctroded,  side,  top  and  corner 
regions,  respectively,  of  the  positive  quadrant  in  act  oruanc  si::: 

LXt 

Figure  J  ot  Ret. lb.  the  e  has  been  suppressed  anc 

3  -  3  cos  ‘i,  3  =3  cos  vo .  3  =3  cos  [i  :  is  co  .  I  . 


The  transcendental  frequency  equations  for  the  trapped  energy  ri¬ 

se  Lucian  may  be  written  in  the  form 

-  -  ,  _S  -  F 

-  tan  a*  =  *  .  v  tan  vb  =  v  ,  o.  .  . 

>  ) 


where 


and 


_S 


a 


,T  = 


3 

“  a 


i.  12) 


and  Che  ei genirequency  a  can  be  iecermined  from  che  dispersion  acua- 
cion  for  Che  eieccrodea  region,  which  cak.es  che  form 

’  2"1  2„.l'  _  •>  _ ' 

c  =  n~~ '  -+n-">  c  -Ml--?-.-.  ■  5 .  ) 

o  n '  n 

ana,  as  noCeti  earLier.  for  cnsuiacmg  films  ^nscaau  >f  eieccrodes 

-  :  1)  -  I  1  ) 

c  is  co  be  replucea  oy  c 

Since  che  main  purpose,  alchough  noc  Che  only  purpose,  of 
chis  secCion  is  for  che  decerminacion  of  che  cemperaCure  dependence 
of  che  effective  piezoe Leccric  ana  iieieccric  conscancs  from  measure- 
mer.es  of  chermaily  compensated  cues,  Lc  is  imp  or  t  an  c  co  noce  chac  in 
addicion  co  u  .  which  appears  in  che  expressions  for  che  crapoea  _ner 
mode,  chere  are  both  u?  and  u^,  which  are  an  order  of  magnicuce 
smaller  chan  u, ,  buc  are  required  in  chis  work  because  che  pure 

L 

Chickness  behavior  is  chermaily  compensaced  for  che  cues  of  primary 
inCeresc  in  chis  seccion.  From  Eqs.(65),  with  i  57)  and  72)  of 
kef. 16,  we  find  chac  co  lowesc  order  in  che  small  wavenumbers  along 
Che  place  che  displacement  field  for  che  S  region  may  be  wriccen  in 
che  form 
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u?-A<l)STl  + 


_s 

•16=  "  C561V;  „  1  .  “"“2  = 

- ITT7 - X:J  SLn  —  e 


(c,  X  -  cetiv)  nnX0  -c  (Xn-i)  ivX, 


“2'[- 

“H' 


_S 


r2=  '  r4LV)  (1)S  nrr  „  ,  ,„(2)S  ^  v  , 

nrr/2h  A+  C°S  2h  X2  L°+  C°S  *2  2h  X2_e 


nrr  ~  *1  (X,  -  *)  i^X- 

htt  1  J 

^  e 


(r5=  '  r3lv)  . (1)S  nr  „  (3)S 

- t~z -  A  .  cos  —  X^  +  lE, 


,S 


ar/2h  A+  C°S  2h  X2  +  lE+  COS  *3  2h  X2je 


nr  ..  1  -=  <V*>  ivX: 


e 

(5.15) 


where 


C<.2)S-i(-l) 


n-1  .  _(2)  ,  x-S  ,  -(2)  .  . 

2  ^r2C  +C12^=  ’  <-r4c  '  c 5 2 ^ 1  ^  ^(1)S 

c  2h)sin(K,nrr/ 2h ) 


n-1  ,  _(3)  ,  .  _S  ,  -(3)  .  N 

(i)c  ?  ^r5c  +C17^=  '  r3c  +c57)lv  ri)s 

Ej_'*E=i(-l)  — 2 rrr — - - - 5 -  A~  ’  (5-lb) 

c  (nr/2h)sin(H^nr/2h) 


and  r,,  r0,  r,  and  are  defined  in  Eqs.(60)  of  Ref. 16  while 
and  are  defined  in  Eqs.(61),  and  (63),,  respectively,  of  Ref. 16. 

It  should  be  noted  that  the  real  part  of  Eqs.(5.15)  is  unders tood ,  and 
we  have  similar  expressions  for  the  ,  T  and  C  regions,  which  clearly 
are  too  cumbersome  to  write  here.  However,  we  note  that  since  we 
have  the  product  of  two  trigonometric  functions  in  the  region,  for 
that  region  the  entire  expressions  for  the  real  part  for  one  direction 
must  be  written  out  explicitly,  thereby  resulting  in  expressions 
twice  as  long. 

The  trapped  energy  solution  presented  in  this  section  is  for 
small  wavenumbers  z  and  v  along  the  plate  and  small  piezoelectric 
coupling.  Consequently,  only  the  thickness  dependence  of  all  elec¬ 
trical  variables  has  been  retained  in  the  treatment^.  In  view  of 
this  and  the  fact  that  the  solution  is  referred  to  the  pLate  axes  and 
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the  eigenvector  triad  of  the  thickness  solution,  the  perturbation 


integral  in  (2.10)  should  be  written  in  the  more  convenient  form 


ll  "  '  I  ’ 


AS<3NvS2,K®v,N  +  2gl,  2  '  A^f,  2^,  2_  dVo  '  (5-17) 


where 


G  .  =  Q,  Q  G'  R  =  3,  R  . 

K3Nv  rv  vo  KvNo  1  lv  v  ' 


(5.18) 


and  1'  is  the  same  as  in  Section  4.  In  Eq.(5.17)  we  have  employed 
the  subscript  1  for  the  dominant  component  of  the  essentially  thick¬ 
ness  trapped  energy  mode  of  interest  in  each  case  because  that  is  in 
accordance  with  the  notation  used  in  Ref. 16.  Consequently,  we  must 
renumber  accordingly  for  each  different  dominant  mode  in  a  given 
plate.  For  the  geometry  of  the  trapped  energy  resonator  considered 
here  from  (5.17)  we  obtain 

h  i  b  = 

H  =  -  4  '  dxJ  r  dX.f  f  (?  +  C)dX,  +  [  (<?T+CT)  duO 

i.  KJ  «jl\0  Jo  J  / 

-h  o  o  b 

fa  => 

(’/(’S2  r  C  r  l 

+  dX-i  1  (9  +C  )  dX  +  !  (9  +  C  )  dX„  .  (5.19) 

J  i\  J  J  o  J_ 

>j  o  b 


wnere 


K3Nv02,K°v,N 


C  =  2£R  f 


-  4Cf 


(5. 20) 


When  the  data  is  available  Eqs.(5.19)  and  (2.9)  are  to  be  used 
to  obtain  4R^  and  LQ  (if  not  known  from  other  measurements)  for 
different  trapped  energy  modes,  including  thickness  overtones,  of 
different  thermally  compensated  cuts.  From  the  AR^  obtained  from 


can 


different  orientations,  (5.18),  and  (4.13),  ^  and  .LR 

readily  be  obtained.  In  a  similar  way  from  the  4C  and  C4.14).  j. 

and  4^2  can  readily  be  obtained.  If  they  are  known  this  can  be  used 

as  a  check.  The  values  of  the  elastic  constants  cT  w  will  be  refined 

2LvMo 

using  the  trapped  energy  analysis  bv  calculating  M  and  P  from  the 

n  n 

previous  determination  of  the  constants  and  using  Eq.i5.I4)  to 
refine  c^  for  the  constant  in  question.  This  will  be  lone  for  ail 
the  modes  to  refine  all  the  constants.  This  will  be  continued 
recursively  until  there  is  no  change  to  a  given  number  cf  significant 
digits.  Then  the  piezoelectric  constants  will  be  refined  from  the 
trapped  energy  analysis  it  necessary.  Finally,  redundant  checks  will 
be  made  on  all  coefficients. 
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